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Abstract 

We obtain the dipolar anisotropies in the arrival directions of ultra-high energy cosmic ray 
nuclei diffusing from nearby extragalactic sources. We consider mixed-composition scenarios 
in which different cosmic ray nuclei are accelerated up to the same maximum rigidity, so that 
E < with Z the atomic number and i?£iax the maximum proton energy. We adopt 

E^^^ ~ 6 EeV so as to account for an increasingly heavier composition above the ankle. 
We obtain the anisotropies through Monte Carlo simulations that implement the cosmic ray 
diffusion in extragalactic turbulent fields as well as the effects of photo-disintegrations and 
other energy losses. Dipolar anisotropies at the level of 5 to 10% at energies ~ 10 EeV are 
predicted for plausible values of the source density and magnetic fields. 


1 Introduction 

In a previous work, hereafter referred to as paper I [Q , we studied the expected dipolar anisotropies 
of ultra-high energy cosmic rays (UHECRs) diffusing from nearby extragalactic sources, under the 
assumption that the cosmic rays were protons. In this work we extend this analysis to the case 
of heavier nuclei, which is of relevance due to the indications provided by measurements with the 
Pierre Auger Observatory of a transition towards a heavier composition at energies larger than a 
few EeV |^, ^ (where I EeV = 10^® eV). 

We will consider scenarios in which a mixed composition is present at the sources and such 
that the acceleration of the different nuclei occurs up to the same maximum rigidity, as expected 
if the acceleration process is of electromagnetic nature. Assuming for simplicity a sharp cutoff, the 
cosmic ray (CR) energies satisfy in this case E < where Z is the atomic number and 

the maximum proton energy adopted. This would then naturally give rise to a transition towards 
heavier elements for energies E > E^^^. We will adopt hereafter E^^^ ~ 6 EeV, the so-called low 
rigidity scenario [^, ||, |^, for which the transition to heavier elements happens just above the ankle 
of the spectrum (which is the observed transition towards a harder spectrum taking place at an 
energy ~ 5 EeV). In this case heavy nuclei such as Fe, which have Z = 26, can reach energies of 
order 150 EeV, comparable to the highest CR energies observed. Note that given the relatively 
small maximum rigidities considered, the conditions required for the properties of the sources are 
less extreme in these scenarios than in the ones in which protons are accelerated up to beyond 100 
EeV. 

An important observable which may help to discriminate between different scenarios is the 
anisotropy in the arrival direction distribution at different angular scales. In particular, the low 
level of anisotropies observed at small angular scales at the highest energies Si {E > 40 EeV) is 
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somewhat unexpected in proton scenarios, for which the deflections in galactic and extragalactic 
magnetic fields would be of only a few degrees. Some hints of intermediate angular scale excesses 
appearing at the highest energies have been reported by the Telescope Array Q (a hot spot on a 
20° radius circle for E > 57 EeV) and by the Auger Observatory Q (an excess on a 12° radius 
circle close to the Super-Galactic Plane and centered at about 18° from the direction of Centaurus 
A, for E > 54 EeV). Finally, searches for anisotropies on large angular scales, which look for 
the presence of a dipole or a quadrupole, are of interest since they could arise in cases in which 
deflections are large and cosmic rays diffuse, and eventually also in the case of quasi-rectilinear 
propagation due to the anisotropic distribution of the nearby sources. The latest results from the 
Pierre Auger Observatory indicate some hints of a ^ 7% dipolar anisotropy for energies E > 8 
EeV @. 


2 Effects of background radiations upon UHECR nuclei 


One of the difficulties that appear when trying to obtain predictions in scenarios involving UHECR 
nuclei is that during their propagation they interact with the radiation backgrounds, both the 
cosmic microwave background (CMB) as well as the extragalactic background light (EBL) which 
contributes mostly at infrared and visible wavelengths. As a result, nuclei can photo-disintegrate 
and change their mass. In the case of protons one just needs to solve an equation describing the 
energy losses, what leads to changes in the Lorentz factor P due to redshift effects as well as the 
photo-pion and pair creation losses. Instead, in the case of nuclei one also needs to consider the 
change in the nuclear mass by the stochastic effects of photo-disintegrations. 

Since we will consider scenarios with relatively small maximum energies, E < 6Z EeV, and 
we will focus on CR energies near and above the ankle, E > 4 EeV, then the relevant Lorentz 
factors for this work will be 7 x 10^ < P < 6 x 10®. In particular, this implies that photo-pion 
production processes will be negligible, since for interactions with the CMB the pion production 
threshold would require values P > 4 x 10^®, while the photo-pion production is negligible in the 
case of interactions with the EBL. 

For the photo-disintegration processes we will adopt the cross section parameterizations given 
by Puget, Stecker and Bredekamp (PSB) |^, with the improved threshold energies reported in 
[|l^ . Given the maximum energies considered, E < 6Z EeV, and the fact that we will only be 
interested in energies E > 4 EeV, we can also safely neglect the secondary nucleons emitted in 
the photo-disintegration process, since they will have energies smaller than E’^^^ZjA ~ 3 EeV. 
This simplifies the analyses since we will just need to follow the leading fragment after the photo¬ 
disintegration. We will also consider that for any given nuclear mass we can assume that the 
relevant nucleus is the one corresponding to the most stable isotope. 

The rate at which a nucleus of mass number A photo-disintegrates with the emission of i 
nucleons is 




1 d£ dn^ 



de' e'a A, 


( 1 ) 


where the cross sections aA,i were parameterized in p4| |. The photon number density dn-,./de 
includes both the GMB and the EBL contributions. An important ingredient in these studies is 
the modelling of the EBL and its redshift evolution, for which we will adopt the values reported 
in the work of Inoue et al. jl^ . 

To describe the average mass loss of a given nucleus it is convenient to introduce the effective 
rate 

^A,ef f — (2) 
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Figure 1: Attenuation lengths for different nuclei (He, C, Si and Fe) and for two values of the 
redshift, z = Q (left) and z = 0.5 (right), as a function of the Lorentz factor F. 


In terms of this rate one has that the average change in the mass number after traversing a distance 
dx is 

dA = -RA,eff dx. (3) 

The attenuation length with respect to the photo-disintegrations can be defined as 


^pd — 


A dx 


A 

RA,eff 


(4) 


We plot in Fig. 1 this attenuation length for different nuclei (He, C, Si and Fe) and for two values 
of the redshift, z = 0 and 0.5, as a function of the Lorentz factor F. Also shown in the plot is the 
attenuation length for pair creation processes, Xpair, computed using the fit obtained in paper I 
for the case of protons and using the scaling relation Ap„j^(r) ~ AAp(jj^(F)/Z^, which holds up to 
small effects associated to Coulomb corrections. 

Using that E ~ FArup, with rup the proton mass, and computing the evolution as a function 
of redshift so as to account more directly for cosmological effects, one has that 



1 dE IdT 1 dA 

(5) 


E dz F dz A dz ’ 

where 

1 dr dt 

T dz ~ ^dz 

■ I I ■ 

_ ^pair ^ad _ 

(6) 

and 

1 dA _ 
A d7 

dt I 
dz Xpd 

(7) 

In the previous expressions 

dt 

I 

(8) 


dz (I-hz)H(z)’ 

with the Hubble constant being 





H{z) — (1 + z)^^m + (9) 
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We will adopt present values Hq = 70 kms“^ Mpc“^ for the Hubble constant, = 0.3 for the 
matter contribution and Ha = 0.7 for the vacuum energy contribution. The attenuation length for 
the adiabatic losses is Xad = c/H(z). Note that the attenuation lengths introduced above have the 
meaning of an energy loss length, i.e. represent the distance over which the energy would change 
to 1/e of its initial value if the attenuation lengths were to remain constant. On the other hand, 
the typical length for the emission of one nucleon by photo-disintegration is Xpd/A. 

In order to account for the stochastic nature of the photo-disintegration processes while follow¬ 
ing the propagation of a given nucleus, we will assume that the probability of emitting one nucleon^ 
in a given step dx is just [1 —exp(— Hdx/Apd)]- In this way we will consider that individual nucleons 
are emitted stochastically with the effective rate that already accounts also for the probability of 
multi-nucleon emission. 


3 Diffusion of UHECR nuclei in turbulent magnetic fields 

At the energies considered in this work, E > A EeV, CRs are most likely of extragalactic origin, and 
our aim here is to study the effects that could arise if they diffuse from their sources in the presence 
of sizable turbulent extragalactic magnetic fields present in the Local Supercluster. Extragalactic 
magnetic fields are expected to be produced either by the evolution of primordial seeds during the 
process of structure formation or as a result of galactic outflows. Their strength should then be 
enhanced in denser regions while suppressed in the voids. A turbulent component in the Local 
Supercluster with RMS values of 1-10 nG, with typical coherence lengths Ic ^ 1 Mpc, is often 
considered, and we will here adopt the illustrative values of R = 1 nG and Ic = ^ Mpc, taken to 
be homogeneous over space. In Section 7 we will analyze the dependence of the anisotropy on the 
assumed strength of the magnetic field. Note that since GRs will be diffusing in this scenario, the 
most relevant magnetic field effects will be those affecting relatively nearby sources, i.e. within 
^ 200 Mpc. Contributions from sources farther away will be suppressed either by a magnetic 
horizon effect at low energies or by the interactions with the radiation backgrounds at the highest 
energies. For definiteness we will consider that the turbulent field has a Kolmogorov spectrum, 
so that the magnetic field energy density scales as u}{k) oc k~^/^ in Fourier space. An important 
quantity (see paper I for details) to characterize the diffusion effects is the critical energy Ec 
corresponding to the value for which the Larmor radius tl becomes equal to Ic, given by 

B I 

Ec = ZeBlc ~ 0.9Z—EeV. (10) 

nG Mpc 

Particles with E < E^ undergo large deflections caused by the magnetic field modes with scales 
comparable to the Larmor radius, giving rise to the so-called resonant diffusion. Instead, particles 
with E > Ec are in the regime in which the deflections over each coherence length, S ~ Ic/i'L, 
are small. Hence, only after the GRs traverse several coherence domains, over a distance scale 
known as the diffusion length Ijj, the overall deflection becomes ^ 1 rad. At distances larger than 
Id = 3R/c spatial diffusion takes place, characterized by the (isotropic) diffusion coefficient D. 
In the regime with E > Ec there is then non-resonant diffusion if the source distance Vs is such 
that Tg ^ Id, while if fs < Id, which is the case at sufficiently high energy, the propagation is 
quasi-rectilinear. 

For the parameters considered in this work we have Ec ~ 0.9Z EeV at present and hence 
for R > 4 EeV only the heavy nuclei may be in the resonant regime. Note that the coherence 
length typically has a redshift dependence lc{z) = lc{9)/{^ + z) while the magnetic field is generally 
enhanced for higher redshifts as B{z) ~ R(0)(1 -I- where the parameter /x was introduced in 

^In the particular case of ^Be we account for the fact that the photo-disintegration produces two alpha particles 
and a nucleon, and for simplicity we assume that the two alpha particles follow similar trajectories afterwards. 
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[p^ to account for MHD effects. Following we will assume /r = 1, in which case the critical 
energy actually results independent of redshift. 

The observed distribution of arrival directions of nuclei that originate from an individual source 
depends on the magnetic field they traverse, the energy of the particles and the distance to the 
source. In the diffusion regime the dipolar anisotropy is given by A = SD'Vn/n, where n is the 
spatial density of particles. In paper I we used this expression to compute A, with n the solution of 
the diffusion equation for protons, including energy losses during propagation, obtained following 
Berezinsky and Gazizov EM We also showed that the dipolar anisotropy can alternatively be 
computed through numerical simulations of the trajectories of particles propagating in a turbulent 
magnetic field in terms of the average of the cosine of the angle 9 between the final CR velocity 
and the vector describing the particle position with respect to the source. The dipole amplitude 
is A = 3(cos0). This method has the advantage that it can be used for all energies, encompassing 
the transition from the diffusive propagation of low rigidity particles, that induces small amplitude 
anisotropies, to the quasi-rectilinear regime at high rigidities, that leads to large anisotropies. We 
implemented this method through a numerical integration of the stochastic differential equation 
that describes the scattering of UHECRs in a turbulent magnetic field 0 MM- 

dui = -^n^cdt + -^=PijdWj, ( 11 ) 

l-D V‘D 

where Pij = {6ij — rurij) is the projection tensor onto the plane orthogonal to h = (ni,n 2 ,n 3 ) 
(the direction of the CR velocity), repeated indices are summed, and (dVFi, dW 2 , dWs) are three 
Wiener processes such that (dWi) = 0 and {dWidWj) = cdtSij. 

In the present paper we use the stochastic method across the entire range of rigidities considered 
to obtain cos 9 for each of the propagated particles. We use as input iov Id = 2>D/c the diffusion 
coefficient D{E) evaluated in paper I through numerical integration of the Lorentz equation for 
the trajectories of charged particles in a homogeneous turbulent magnetic field. The analytic fit 
to the results for the case of Kolmogorov turbulence (that were comparable to those in ref. jl^ ) is 
given by 


D{E) = 




( 12 ) 


If energy losses and cosmological expansion effects are neglected, the dipole anisotropy evaluated 
as the average A = 3(cos0) over a large set of trajectories obtained with a numerical integration of 
equation 0 coincides with A = SD{E)/crs for Vg > Id, as expected in the diffusion regime. In 
the case of much smaller than Id it coincides with the approximation A/3 ~ 1 — crs/{9D{E)) 
valid for the quasi-rectilinear regime (an approximate expression for the anisotropy valid for all 
source distances and energies is given by eq. (21) in paper I). The use of this method allows to 
evaluate the anisotropy across the different regimes, accounting also for energy losses and nuclear 
fragmentation, as we do in the next sections. 


4 Results for individual sources and for each type of nuclei 

We first compute the dipolar anisotropy from an individual source that injects just one type of 
nucleus, as a function of the source distance and the arrival energy. To cover the relevant range of 
possible masses we consider five different representative nuclei: p. He, C, Si and Fe. In each case 
we propagate particles with energies ranging from Emin = 4 EeV up to Emax = dZ EeV, with a 
spectral distribution E~'^ and a sharp cutoff at Emax- We inject at the source a constant number 
of particles per unit of time, neglecting for simplicity source evolution effects, and we integrate 
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Figure 2: Left: dipole amplitude as a function of the arrival energy for a source at distance of 
25 Mpc (circles), 150 Mpc (squares) and 300 Mpc (plus signs) emitting He, C or Fe nuclei. Right: 
mean mass as a function of the final energy for the same cases. A magnetic field strength B = 1 nG 
is adopted. 


the contributions from a maximum redshift Zmax = 1 down to the present, z = 0 (note that in 
the diffusion regime the redshift has the meaning of time rather than distance). We follow the 
evolution of the mass number (and charge) and of the Lorentz factor of the particle as a function of 
the redshift as described in Section 2. We also follow the spatial trajectory as described in Section 
3, solving equation dO) numerically, and evaluating cos 9 for each propagated particle. When the 
particle arrives to z = 0 we record its energy, mass, distance from the source and the angle between 
the source direction and the CR arrival direction. We can then obtain the dipolar amplitude of 
nuclei by computing (cos 9) averaged over the set of particles with specified arrival energy and 
source distance. 

For example, we show in the left panel of Figure || the dipole amplitude for injected He, C and 
Fe nuclei coming from sources at three different distances, 25 Mpc, 150 Mpc and 300 Mpc, as a 
function of the arrival energy. Comparing the different curves we see that for a given energy the 
lighter particles lead to larger anisotropies. This arises because having larger rigidities they suffer 
smaller deflections than the heavier ones. For the closer distance displayed of 25 Mpc, for which 
the energy losses are less relevant, the anisotropies of the different components are quite similar 
for the same rigidity. For instance, shifting the Fe curve to the left by rescaling the energy by a 
factor 6/26 would closely match the results obtained for C. 

When considering sources farther away the anisotropies first diminish with increasing distance, 
being smaller e.g. for 150 Mpc than for 25 Mpc, since the diffusive flux becomes more isotropic 
farther away from the source, but then the anisotropy may increase again (see for instance the 
300 Mpc curve for Fe). This is the consequence of two effects. On one side the magnetic horizon^ 
enhances the gradient of the CR density at very large distances, and on the other hand, due 
to photo-disintegrations the contribution from lighter particles becomes significant for faraway 
sources. Since the photo-disintegrations reduce the energy of the leading fragment (by reducing its 
mass) but they leave essentially unchanged the particle rigidity (since A/Z ~ 2 for nuclei heavier 

^The magnetic horizon is essentially the distance from a source that a CR can travel diffusively in a time 
comparable to the age of the Universe. Hence, the density of CRs produced by the source should be strongly 
suppressed beyond the (energy dependent) magnetic horizon distance. 
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Figure 3: Density of particles times (in arbitrary units) as a function of the distance to the 
source. The left panel is for C injected nuclei for energies of 4, 8, 17 and 34 EeV (from top to 
bottom). The right panel is for injected Fe nuclei and for energies of 4, 8, 17, 34 and 70 EeV. 


than H), this implies that at a given energy the average rigidity increases farther from the source, 
and this can lead to an increase in the resulting anisotropies with distance. In particular, a peculiar 
feature that is observed for intermediate nuclei, as in the C case, is a kink in the dipole amplitude 
curves appearing at ^ 10 EeV. This can be understood by looking at the right panel of Figure |[ 
where the mean mass of the arriving particles is plotted for the same injected nuclei and source 
distances considered in the left panel. The sudden drop in the mean mass for the injected C 
case appearing for energies below ^ 10 EeV is due to the effects of the Be disintegration to two 
alpha particles plus a nucleon. Since the maximum C energy considered is 36 EeV, in this case 
the maximum energy for the alpha particle resulting from the photo-dissociation process would be 
about 12 EeV, and this is why the average mass number has a significant drop below this energy. 
Since the alpha particles appearing below 12 EeV have a relatively high rigidity, their contribution 
to the anisotropies is enhanced, explaining the rise in anisotropy observed below 12 EeV, and we 
also see that in this energy range the anisotropy from the faraway source would be comparable to 
that obtained in the case of injected He nuclei. 

From the results of the simulations it is also possible to obtain the density of particles as a 
function of the distance to the source for different values of the energy. We show in Figure ^ the 
results for injection of C and Fe nuclei. At the highest energies, the density decreases as n oc r“^, 
as expected from quasi-rectilinear propagation, but particles quickly loose energy and cannot arrive 
from far away sources (farther than ~ 50 Mpc for 34 EeV carbons or ~ 150 Mpc for 70 EeV irons). 
For smaller energies the losses are less important and particles can arrive from sources farther 
away. Moreover, an enhancement of the density with respect to the law is evident at short 
distances, characteristic of the diffusive propagation, followed then by a steep decrease at large 
distances due to the magnetic horizon effect. 

5 Large scale anisotropies from a distribution of sources 

In the previous section the dipolar anisotropy from individual sources in the presence of a turbulent 
magnetic field was computed as a function of the source distance, type of nucleus injected and of 
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Figure 4: Mean total dipole amplitude and its dispersion as a function of the energy for different 
injected nuclei and a density of sources p = 10“® Mpc“^ (red solid lines) and p = 10“"*^ Mpc“^ 
(black dashed lines). 


the CR energy. In a realistic situation the total cosmic ray flux will probably originate from a set 
of several (or many) sources, and probably accelerating a mixture of different nuclei. The total 
dipolar component of the flux will mainly depend on the location and intensities of the nearest 
sources and on whether there is an inhomogeneous distribution of the sources at large scales. If 
there are several sources contributing to the flux, the dipolar anisotropy can be obtained from the 
superposition of the individual source dipoles. 

We first consider in this Section the case in which just one type of nucleus, j, is injected by all 
the sources. The corresponding dipolar anisotropy is 

^ U) 

A^^\E) = J2—iE)^^\E), (13) 

where N is the number of sources giving a non-negligible contribution to the flux at energy E, 
nf'^/nt measures the fraction of the total flux coming from the f-th source and aI^\e) is the 
dipole anisotropy of the flux from source i as computed in the previous section for a given nucleus 

j- 

In order to estimate the relative contribution to the anisotropy of the different sources as a 
function of their distance to the observer we will make the simplifying assumption that the sources 
are steady and have equal intrinsic intensities, so that for each energy the relative contribution to 
the flux from different sources will only depend on the distance to the source r^. 

The fact that the sources are distributed in different sky directions means that the vector sum 
in eq. generally leads to a smaller dipole amplitude when many sources contribute. In the 
case that only few sources are relevant, the direction of these particular sources will determine the 
dipolar anisotropy, while if many sources are relevant, the overall large scale distribution of the 
sources, in particular whether the distribution has a non-vanishing dipole component, can have a 
significant effect. 

In order to quantify the total amplitude of the dipolar anisotropy we perform some simple 
simulations considering first a homogeneous distribution of sources. Starting with one source at a 
random direction in the sky, that represents the closest source, we subsequently add new sources 
in random directions and compute the new total dipolar anisotropy using eq. di)- The radial 
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distances from the observer to the sources are taken as the mean expected value for the i-th closest 
source in an homogeneous distribution, that is given by (r^) = (3/47rp)^/^r(i + l/3)/(i — 1)!, where 
p is the density of sources. We show in Figure ^ the mean amplitude of the dipole obtained in 1000 
simulations for two different values of the source density: p = 10“® Mpc“^, for which the closest 
source is at a mean distance (ri) ~ 25 Mpc (solid lines), and p = 10“^ Mpc“^, for which (ri) ~ 11 
Mpc (dashed lines) and for the five nuclear species considered. For a larger density of sources the 
flux from each source is smaller, leading to smaller anisotropies. 


6 Dipolar anisotropy in a mixed-composition scenario 


We can now consider a more realistic scenario in which a distribution of sources, that we again 
assume to have equal intensity and to be homogeneously distributed, accelerates cosmic rays with 
a mixed composition. For simplicity we will consider the same composition for all the sources, 
that we parameterize by the fraction fj of the differential energy flux contributed by nucleus j at 
low energies (below the proton spectrum cutoff) of each of the five representative nuclei that we 
considered in the previous section. The expected dipolar anisotropy is then obtained by a weighted 
superposition of the amplitudes corresponding to each nuclei, 

A{E)=Y,f,^A\^\E), (14) 

. . 


where nt{E) = ^ ^ (E) and (E) is the dipolar anisotropy of source i assuming that only 

the individual nuclei j contribute, as obtained in the previous section. 

As an example we will show the results for a particular set of nuclei fractions at the sources that 
lead to a cosmic ray spectrum similar to that measured by the Pierre Auger Observatory 120 . The 
left panel of Figure ^ shows the overall spectrum of particles reaching the observer from sources 
accelerating particles with a spectrum E~^ up to maximum energy E^ax = bZ EeV and injected 
fractions fp = fne = fsi = 0-19,/c = 0.4 and fpe = 0.03, propagating through a homogeneous 
turbulent magnetic field with Kolmogorov spectrum, coherence length of 1 Mpc and root mean 
square amplitude of 1 nG. Given the assumed E~^ spectrum and the sharp cutoffs adopted the 
spectrum obtained with the chosen fractions is indeed quite acceptable. Also shown in Figure 
are the separate contribution to the spectrum from the three ranges of masses, corresponding to 
A < 4, 5 < A < 16 and 17 < A < 56. 

The composition as a function of the energy also follows the general trends measured by the 
Pierre Auger Observatory, as it can be seen in the right panel of Figure |^, where the mean (In A) 
and variance V (In A) are shown for the model and are compared to the values obtained from the 
data adopting the EPOS-LHC hadronic interaction model. 

The total dipolar anisotropy from homogeneously distributed sources with spatial density p = 
10“® Mpc“^ and p — 10“"^ Mpc“^ and accelerating the nuclei with the above mentioned fractions 
is shown in the left panel of Figure ^ as a function of the energy. Here the mean and dispersion of 
the dipolar amplitude obtained from 500 simulations are plotted. The sharp features appearing in 
the curves at the energies corresponding to maximum acceleration for the five elements considered 
are an artifact of the sharp cutoffs adopted, and are expected to be smoothed out in scenarios with 
softer suppressions, such as exponential ones. 

In the right panel of Figure ^ we plot the dipole amplitude separately for different mass groups. 
We divide the composition observed at Earth in three sets, A < 4, 5 < A < 16 and 17 < A < 56, 
the same ranges used to illustrate the different contributions to the spectrum in the left panel of 
Figure |[ We show the dipolar anisotropy amplitude for each mass group as a function of the 
energy for the case of a homogeneous distribution of sources with density p = 10“^ Mpc“^. The 
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Figure 5: Left: Spectrum of CRs. The data points are from the Pierre Auger Observatory |^, the 
black solid line is for a homogeneous distribution of sources accelerating a mixed nuclei composition 
with ifp, Jhe, fc, fsiT fpe) = (0.19,0.19,0.4,0.19,0.03). The red dashed line is the contribution 
from nuclei with A < 4, the blue dotted line is that from nuclei with 5 < A < 16 and the purple 
dot-dashed line is for 17 < A < 56. Right: composition at Earth in terms of (InA) and its variance, 
R(lnA). The data points are from the Pierre Auger Observatory |Q based on an interpretation of 
air-shower data using EPOS-LHC. Error bars correspond to the statistical uncertainties and color 
bands to the systematic ones. The solid lines are the predictions for the same model used in the 
left panel. 


amplitude in the lighter components is enhanced relative to the total one since the heavier fraction 
is more isotropically distributed at lower energies. 

The previous results hold for the case of homogeneously distributed sources. In the case that 
the sources themselves have an inhomogeneous distribution around the observer, in particular if 
their distribution has a non-vanishing dipole, a further contribution to the anisotropy is expected. 
For instance, if the distribution of nearby cosmic ray sources follows the local distribution of matter 
a non-vanishing dipole is expected. The dipolar component of the matter distribution is indeed 
known to be responsible for the Local Group peculiar velocity with respect to the rest frame of 
the CMB, that actually gives rise to the observed CMB dipole. The dipolar component of the 
mass distribution in our neighborhood has been estimated using different catalogs of galaxies, 
as for example the 2 Micron All-Sky Redshift Survey (2MRS), showing that the resulting dipole 
seemingly converges when sources up to a distance ~ 90 Mpc are included |2^ . 

The effect of the local inhomogeneity of the source distribution in the predicted large scale 
anisotropies can be included in the simulations by choosing the positions of the sources in our 
neighborhood from some catalog representing the local distribution of matter. To describe the local 
distribution of matter we use a volume limited subsample^ of the 2MRS catalog up to 100 Mpc 
. We have then selected the position of the required number of sources (according to the density 
considered) from this subsample of 2MRS galaxies. On the other hand, the locations of sources 
farther away were assumed to be isotropically distributed. We show in Figure ^ the change in the 
mean dipole amplitude when the inhomogeneous source distribution is considered, adopting the 
same mixed-composition scenario analyzed above, a density p = 10“^ Mpc“^ and a turbulent held 
of R = 1 nG. Compared to the homogeneous case, an enhancement of the dipole amplitude by a 

^Considering only objects with d < 100 Mpc and absolute magnitude in the K band Mjy < —23.4. 
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Figure 6: Left: Mean total dipole amplitude and its dispersion as a function of the energy for a 
homogeneous distribution of sources with density p = 10“^ Mpc“^ (red solid lines) and p = 10“^ 
Mpc“^ (black dashed lines) accelerating a mixed nuclei composition with (/p, Jme, fc, fsi, fpe) = 
(0.19,0.19,0.4,0.19,0.03). Right: mean dipole amplitude discriminating by composition at Earth, 
in the mass ranges as indicated by the labels, for the case of sources with density p = 10“"* Mpc“^. 


factor of about 2 is observed. 


7 Anisotropies and the strength of the magnetic field 

In the previous sections we computed the resulting anisotropies, both in the case of a single source 
as well as for a superposition of sources, adopting a reference value of 1 uG for the RMS turbulent 
extragalactic magnetic field B. We now discuss how the predictions would be affected for different 
assumed strengths of the magnetic field. 

The left panel of fig. ^ shows the energy dependence of the anisotropy of an individual source 
of Fe nuclei at a distance of 25 Mpc, under the assumptions oi B — 0.1, 1 and 3nG. As expected, 
the resulting anisotropy values in the different cases are comparable if one rescales the energies by 
the inverse of the assumed magnetic field strength. This results because the diffusion coefficient D 
is just a function of E /and the critical energy is proportional to B. To better appreciate this 
behavior we show in the right panel the resulting anisotropies as a function oi EjB for sources 
at 25 and 100 Mpc. The approximate scaling is apparent, but is however violated by the effects 
of interactions with the background radiation, since the interaction cross sections just depend 
on E. Also the fact that the surviving leading fragment can have a slightly different rigidity 
(since AjZ may vary) can affect the final anisotropy. We see then that the overall impact of 
photo-disintegrations is more pronounced for sources farther away and the anisotropies for larger 
magnetic fields are larger than what would be suggested by the simple scaling. This increase is 
mostly due to the fact that for stronger B fields the distance travelled by the GRs from the source 
is larger, and hence the photo-disintegration effects are stronger, so that for a given arrival energy 
the actual rigidity of the arriving particles will be larger, and hence the anisotropy will be stronger. 

Turning now to the case of a superposition of sources, here an interesting effect appears that 
leads to a reduced dependence of the dipolar amplitude with respect to the strength of the magnetic 
field. Indeed, for a given nearby source (closer than the magnetic horizon), as the value of B is 
increased the individual contribution of that source to the GR density increases because its density 
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Figure 7: Mean total dipole amplitude and its dispersion as a function of the energy for a distri¬ 
bution of sources with density p = 10“^ Mpc“^ homogeneously distributed in space (black dashed 
line) and following the local matter distribution (blue solid line) that accelerate a mixed nuclei 
composition with the same injected fractions as in Figure 


is enhanced by the diffusion while that from far away sources is suppressed. On the other hand, 
the value of the dipolar component of its anisotropy decreases in such a way that both changes 
compensate each other to a large extent. In particular, in the diffusion regime and neglecting 
interaction and redshift effects, one has that the CR density from one source satisfies 


Q 

* dTTriZ?’ 

with D the diffusion coefficient, while its contribution to the dipole is 


(15) 


A, = 


ZD 

Ti 


We see then that for a superposition of sources 


^ ntot 



ZQfj 
47rr| ’ 


(16) 


(17) 


which is independent of the value of D, and thus of the magnitude of B, and the result is similar 
to what would be obtained in the absence of magnetic fields. We have here used that, due to the 
propagation theorem | |2^ , the total contribution to the CR density ntot is independent of the value 
of B, which is valid as long as the source density is uniform and the inter-source distance is smaller 
than the diffusion and interaction lengths. We have checked numerically that the anisotropy result 
in eq. also holds, as long as interaction and redshift effects are ignored, in the regime of 

quasi-rectilinear propagation. 
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Figure 8: Left: Dipolar anisotropy as a function of the arrival energy for an Fe source at 25 Mpc 
in the presence of a turbulent magnetic field of 0.1 nG, 1 nG and 3 nG. Right: In blue (small 
symbols) the same results as in the left plot but displayed as a function of E/B. In red (large 
symbols) the corresponding results for a source at 100 Mpc. 



E [EeV] 


Figure 9: Mean total dipole amplitude as a function of the energy for a homogeneous distri¬ 
bution of sources with density p = 10“® Mpc“^ accelerating a mixed nuclei composition with 
(/pj Ihe, fc, fsi, fpe) = (0.19, 0.19, 0.4,0.19,0.03) in the presence of a turbulent magnetic field of 
0.1 nG, 1 nG and 3 nG. 

In fig. ^ we show the resulting anisotropies, including interaction losses and redshift effects, 
for the case of a uniform source distribution with density 10“^ Mpc“^, for the same composition 
mixture considered in fig. || and for the three values B — 0.1, 1 and 3nG. We see that differences 
in the dipolar component of the anisotropy by up to a factor of 2 can arise due to the interaction 
effects. 
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It is important to note however that when the deflections from the nearby sources become 
small, i.e. at high rigidities, the CRs arriving from those sources will be quite localized in the 
sky. Hence, although the dipolar component of the arrival direction distribution may have only a 
moderate dependence on the actual strength of the B field, the detailed distribution will be quite 
different in both cases. In particular, higher harmonics will become more sizeable for weaker B 
fields and the individual contributions from the nearby sources would become distinguishable in 
this case. 

Finally, let us mention that the additional deflections of the CRs due to the galactic magnetic 
held, not included in this work, would modify the CR dipole amplitude and direction as well as 
generate higher order multipoles in the arrival directions distribution but for the energies 
considered here the induced change in the dipole amplitude is expected to be small. 

8 Conclusions 

The low level of anisotropies in the distribution of arrival directions of UHECRs ill challenges 
scenarios with a relatively small density of cosmic ray sources and a light composition. In addition, 
there are hints of a dipole anisotropy with amplitude around 7% for energies above 8 EeV 0. 

We have shown in paper I |Q, where we assumed the CRs to be protons, that a dipole anisotropy 
with amplitude between 3% and 10% at 10 EeV is expected in scenarios in which CRs that origi¬ 
nate in extragalactic sources diffuse in turbulent magnetic fields. This result was derived for source 
densities between 10“'^ Mpc“^ and 10“^ Mpc“^, magnetic fields with strength of order 1 nG and 
coherence length of the order of 1 Mpc. Here we extended the calculation of the predicted dipole 
anisotropy as a function of CR energy to the case of nuclei. We developed a method to evaluate 
the anisotropy across different regimes, from spatial diffusion to quasi-rectilinear propagation, ac¬ 
counting for energy losses and nuclei fragmentation. We applied the method to mixed-composition 
scenarios in which different cosmic ray nuclei are accelerated up to the same maximum rigidity. 
With sufficient magnetic turbulence, each component contributes significantly only in the energy 
range between its low-energy cut-off due to a magnetic horizon effect and its maximum acceleration 
limit. We illustrated the method in a scenario with maximum acceleration energy Emax = 6Z EeV, 
and fractions of different nuclei fp = fue = fsi = 0.19, fc = 0.4 and fpe = 0.03 injected with an 
E~‘^ spectrum, which provides an acceptable match to the observed spectrum [0| and composition 
1^ derived from measurements by the Pierre Auger Observatory. 

Naively one could expect a smaller dipole anisotropy at a given energy in a scenario with nuclei 
instead of just protons, due to the larger magnetic deflections. We found instead that the mixed- 
composition scenario considered here also leads to dipole amplitudes close to 5-10% at 10 EeV, 
and larger values at higher energy, for the same source densities and magnetic fields as quoted 
above for the proton-only model, and showed that the dependence with the actual value of the 
magnetic field strength is moderate. As can be seen in Figure the dipole amplitude in the arrival 
directions for different injected nuclei are comparable for similar rigidities, even though energy loss 
and fragmentation processes introduce important differences. In the mixed-composition scenario 
considered here, the anisotropy due to a given type of nucleus below the respective Emax is on one 
side enhanced due to the suppression of faraway sources for energies close to E^ax, but on the other 
hand the overall anisotropy is also reduced by the effects of heavier, more isotropic components, 
leading to the trend depicted in Figure ^ The inhomogeneity in the local source distribution 
enhances the predicted dipole amplitude by a factor of about 2 compared to the homogeneous 
case, as shown in Figure |^. Since this mixed-composition model predicts dipole amplitudes in 
the lighter components larger than the total one (at energies below the respective Emax), an 
observatory with capability for mass-discrimination could provide additional tests of this scenario 
and improve the determination of its parameters. 
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